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The redshifted 21-cm background is expected to be a powerful probe of the early Universe, carrying 
both cosmological and astrophysical information from a wide range of redshifts. In particular, the 
power spectrum of fluctuations in the 21-cm brightness temperature is anisotropic due to the line- 
of-sight velocity gradient, which in principle allows for a simple extraction of this information in 
the limit of linear fluctuations. However, recent numerical studies suggest that the 21-cm signal 
is actually rather complex, and its analysis likely depends on detailed model fitting. We present 
the first realistic simulation of the anisotropic 21-cm power spectrum over a wide period of early 
cosmic history. We show that on observable scales, the anisotropy is large and thus measurable at 
most redshifts, and its form tracks the evolution of 21-cm fluctuations as they are produced early 
on by Lyman-a radiation from stars, then switch to X-ray radiation from early heating sources, and 
finally to ionizing radiation from stars. In particular, we predict a redshift window during cosmic 
heating (at z ~ 15), when the anisotropy is small, during which the shape of the 21-cm power 
spectrum on large scales is determined directly by the average radial distribution of the flux from 
X-ray sources. This makes possible a model-independent reconstruction of the X-ray spectrum of 
the earliest sources of cosmic heating. 
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Introduction. The high-redshift intergalactic 
medium (IGM) can be probed by measuring the in¬ 
tensity of photons in the redshifted 21-cm line of the 
hydrogen atom [l!-0|. The 21-cm intensity (expressed 
as a brightness temperature) contains a mixture of 
astrophysical and cosmological information Q from the 
era when the Universe was mostly neutral, i.e., before 
and during the Epoch of Reionization (redshift z > 7). 
Because it is coupled to numerous inhomogeneous 
inputs, the 21-cm signal fluctuates on all scales. Among 
the most important sources of perturbations are those in 
the matter density, 6, the peculiar velocity, and radiative 
backgrounds produced by stars and their remnants; the 
latter include X-rays which heat the IGM, Ly-a photons 
which couple the 21-cm line to the gas temperature, 
and ionizing photons. The various fluctuation sources 
dominate at different epochs throughout cosmic history, 
leading to a succession of peaks in the three-dimensional 
power spectrum of the 21-cm brightness temperature 
fluctuations (hereafter the ”21-cm power spectrum”) 
averaged over the line of sight @43- However, the exact 
size and scale-dependence of the 21-cnr fluctuations 
are still highly uncertain, since they depend on the 
properties of high-redshift astrophysical sources, which 
are poorly constrained due to the lack of observations at 
present. 

Of particular interest is the nature of the first heat¬ 
ing sources which raised the gas temperature above the 
temperature of the cosmic microwave background (CMB) 
radiation. The candidates for these sources are varied 


and include X-ray binaries, mini-quasars, and hot gas in 
the first galaxies, plus more exotic possibilities such as 
decaying dark matter. The spectral energy distribution 
(SED) of X-ray photons emitted in each case can be very 
different, ranging from a soft power-law spectrum js] to 
a hard spectrum that peaks at energies of several keV as 
in the case of X-ray binaries iEi Because the heating 
is inhomogeneous, the character of the SED is imprinted 
in the 21-cm power spectrum, which can be used to con¬ 
strain the average energy of the X-ray photons which heat 
the IGM ]J. j_2 . 


Since 21-cm fluctuations arise from a superposition of 
several different sources of fluctuations, the signal is com¬ 
plex and its interpretation requires detailed modeling. 
Finding additional probes of the effects of various as¬ 
trophysical objects would be a major step forward, es¬ 
pecially more direct, mo del-independent ways to reveal 
their properties. One possible tool is the directional de¬ 
pendence of the 21-cm fluctuations. 


Anisotropy of the power spectrum. The 21-cm 
power spectrum is predicted to be anisotropic due to the 
radial component of the peculiar velocity gradient cre¬ 
ated by structure formation flij - ITii] ] . Specifically, in a re¬ 
gion with a lower velocity gradient than average, a larger 
stretch of atoms along the line of sight contributes to the 
optical depth of the 21-cm line. Within linear theory, 
this effect results in a power spectrum that is a simple 
sum of terms [15] , 


P(k,z) = P ll o(k,z) + 3P IJ ,2(k,z)nl + 5P fJ ,4(k,z)nt , (1) 
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where /ik = cos 0 k in terms of the angle 6 k between the 
wavevector k of a given Fourier mode and the line of sight 


The /independence of the 21-crn power spectrum can 
in principle be used to differentiate between the cosmo¬ 
logical and astrophysical contributions to the fluctuations 
|15j . In particular, since the velocity gradient is directly 
related to density through the continuity equation, P^i 
is proportional to the primordial density power spectrum 
Pg , P^o is the sum of all the isotropic 21-cm fluctua¬ 
tion sources (i.e., those other than the velocity gradient), 
and P ^2 is proportional to the correlation of density with 
the isotropic 21-cm fluctuation sources. However, 21-cm 
fluctuations are non-linear, and recent numerical inves¬ 
tigations during cosmic reionization 17 - 20| suggest that 
it will be difficult to make this separation of terms and 
recover cosmological information from the 21-cm power 
spectrum. 

These investigations focused on the P^i term, which 
carries the primordial cosmological information. This 
term is typically the smallest, and is most easily contam¬ 
inated by non-linearity 21]. In what follows, we rekindle 
the importance and usefulness of the linear-theory poly¬ 
nomial decomposition of the power spectrum (eq. [D by 
focusing on the P ^2 term instead. 

Realistic mock data. While upcoming 21-cm ob¬ 
servations will be hard-pressed to beat other probes in 
terms of basic cosmology, they are likely to reveal to 
us several chapters in the history of stars and galaxies 
for the first time. We use simulated 21-cm observations 


for the first realistic study [22] of the anisotropy in the 
21-cnr power spectrum during the entire period of cos¬ 
mic history relevant for upcoming experiments, from the 
formation of the first substantial population of stars (at 
z ~ 30) until the end of reionization (z ~ 7). We create 
mock data using a hybrid, semi-numerical simulation in 
which the non-linear evolution of the redshifted 21-cm 
signal is followed accounting for its dependence on <5, the 
line of sight velocity gradient, and inhomogeneous X-ray, 
Ly-ct and ionizing radiative backgrounds (e.g., .11] and 
references therein). The simulated volume is 384 3 Mpc 3 
with a 3 Mpc resolution, and assumes that stars form 
via atomic cooling. In addition, we average our results 
over twenty randomly generated sets of initial conditions 
to decrease the statistical errors on large scales. On the 
resolved scales, perturbations in the density and velocity 
gradient fields are linear, but those in the galaxy density 
(and thus in the radiation and 21-crn fields) are enhanced 
(“biased”) and thus non-linear i23j. 

Results and discussion. Among the various terms 
in eq. CD. P^o dominates the usually discussed total 
(angle-averaged) 21-cnr power spectrum, so we focus on 
P ^2 as the most promising additional source of observable 
information. We find that 4 is more difficult both to 
measure and to interpret (so we mostly leave it for future 
work). Now, as noted above, in linear theory, Pjfi would 


be proportional to the cross-correlation of density with 
all the isotropic 21-cm fluctuation sources [24[. While 
the total 21-cnr power spectrum does track the changing 
sources of fluctuations indirectly through changes of mag¬ 
nitude and slope, P” 1 does this much more unambigu¬ 
ously and explicitly through changes of sign [H] . This is 
due to the fact that some sources are positively correlated 
with density, and others are negatively correlated. 

• Tracking cosmic history. As shown in Figure [T] the 
actual observed P ^2 (obtained from fitting the form of 
eq. © to mock data) mostly tracks P^l ', as (going from 
early to late times) it is positive during the epoch of 
Lyman-a fluctuations, goes negative when X-ray heating 
fluctuations take over, then positive again as heating fluc¬ 
tuations continue but the mean IGM temperature rises 
above that of the CMB (a key cosmic milestone termed 
the “heating transition”), and finally back to negative 
once ionization fluctuations come to dominate Eij . 

We also consider a simpler, more robust measure of 
anisotropy that does not require fitting the power spec¬ 
trum to a particular form. We define the ’’anisotropy 
ratio” 


'MvAj Z ) — / T>( 1 \ \ 1 ’ 

(P(k, z) I/Xfclco.5) 


( 2 ) 


which indicates how much P(k, z) increases or decreases 
with /ifc, on average 27]. In particular, when the power 
spectrum has a weak angular dependence, is close to 
zero; on the other hand, when the power spectrum as 
a function of fik has a considerable positive or negative 
tilt, the anisotropy ratio is far from zero and its sign is 
positive or negative, respectively. 

The first conclusion from considering the anisotropy 
ratio (which is also shown in Figure [T]) is that the 
anisotropy is generally large (i.e., |r M | is of order unity) 
and thus potentially observable. (We additionally find in 
our simulation that its measurement is not very sensitive 
to the lowest values of |/ifc|, which are highly foreground 
contaminated 28 -3lj]). More importantly, is rich in 
information. Thanks to the characteristic structure of 
peaks and plateaus, the timing of cosmic events can ap¬ 
proximately be read off the anisotropy ratio. 

Indeed, the sign and general shape of ry, mostly track 
those of P^ 2, since the P„2 term usually dominates the 
angle-dependent terms. There are interesting differences, 
though. When P fl 2 passes through zero (as it is switch¬ 
ing sign), this usually indicates that two different fluctua¬ 
tion sources (one positively and one negatively correlated 
with density fluctuations) are superposing and canceling 
out (at a particular wavenumber), so P^o is also at that 
moment small compared to the velocity gradient term, 
and has a significantly positive value (P)T is always 
positive, and the non-linear r^ reflects this). However, 
when P ^2 = 0 due to the heating transition, ~ 0 as 
well since the heating fluctuations continue to dominate 
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Hard SED, k = 0.2 [Mpc -1 ] 

V [MHz] 

178.6 119.0 71.4 47.6 



Hard SED, k = 0.5 [Mpcf 1 ] 
v [MHz] 

178.6 119.0 71.4 47.6 



Soft SED, k = 0.2 [Mpcf 1 ] 
v [MHz] 

178.6 119.0 71.4 47.6 



Soft SED, k = 0.5 [Mpc -1 ] 
v [MHz] 

178.6 119.0 71.4 47.6 



FIG. 1. The anisotropy ratio r M (top plot within each panel) and coefficient P^{k,z) (bottom plot within each panel). For 
the latter, we compare the actual observable 2 (red) to the theoretical result based on linear separation (black). The 
observed 2 is shown only where the fit (to eq. [lj that produced it works well (i.e., R 2 > 0.9). Positive values of P M 2 or 
Pkf are shown with solid lines, and negative with dotted lines. The results are shown for the cases of a hard X-ray SED (left 
panels) or a soft SED (right panels) at two wavenumbers, k = 0.2 Mpc -1 (top panels) and k = 0.5 Mpc -1 (bottom panels). A 
thin horizontal line in each panel indicates r^ = 0. The thin vertical lines in each panel mark global milestones in the history 
of the Universe; from left to right: the mid-point of cosmic reionization (i.e., when half the cosmic gas has been reionized), the 
redshift of the heating transition, and the peak of the Ly-a coupling era (specifically, the redshift of the peak fluctuation level 
at k = 0.3 during this era). 


through this transition and are not canceled out by a 
second source. Note that the heating transition can be 
recognized (at least at k = 0.2 Mpc^ 1 , for which the 
heating fluctuations clearly dominate around that time) 
as the point when = 0 as it rises with time from neg¬ 
ative to positive values. 


While current observations point to a most likely sce¬ 
nario of heating by X-ray binaries with a hard spectrum 
num, there remains a great uncertainty about X-ray 


sources at such high redshifts. To illustrate how the spec¬ 
trum of the X-ray sources can be probed, we also show re¬ 
sults for a soft power-law SED Q normalized to the same 
total emitted energy. We focus on relatively large (i.e., 
cosmological) scales, and show two different wavenum¬ 
bers in order to illustrate the scale dependence. The 
observable quantities (P M 2 and r^) reflect the fact that 
for the hard SED (compared to the soft one), the heat¬ 


ing transition occurs later, and the heating fluctuations 
are smaller (especially on small scales), so they dominate 
over other sources for a shorter time. The heating is late 
since only a fraction of the energy of the hard X-rays 
is absorbed, and the heating fluctuations are small since 
the hard X-rays mostly come from large distances 111. 


• Measuring the spectrum of early X-ray sources. Rec¬ 
ognizing our ability to track cosmic history using 2 
and/or makes possible an especially promising direct 
measurement. Among the three cosmic periods consid¬ 
ered above, the X-ray heating era is the best target, 
since its source properties are highly uncertain (unlike 
the Lyman-a spectrum, which is predicted to be rela¬ 
tively similar for modern and primordial stars i) , and 
it is also relatively easy to interpret (unlike the bub¬ 
ble structure of reionization that is a result of the small 
mean-free-path of ionizing photons in the neutral IGM). 
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Also, during the era when heating fluctuations dominated 
the 21-cm fluctuations there was a unique moment (the 
heating transition) when the global mean 21-cnr inten¬ 
sity was zero (relative to the CMB). At this time, all the 
other fluctuation sources are suppressed (e.g., to linear 
order they are nullified, since each fluctuation gets mul¬ 
tiplied by the mean 21-cm intensity). This suppression 
includes the velocity gradient effect, so the power spec¬ 
trum is nearly isotropic [25] (i.e., P M 2 and r^ are both 
close to zero). 

The analyses of 21-cm fluctuations due to fluctuations 
in Lyman-a or X-ray radiation are quite similar [B|, [s]. 
When fluctuations are dominated by one of them, in lin¬ 
ear theory the predicted power spectrum can be written 
as 

P(k,*) «i>(M) = [T b (z)Wjr(k,z)] 2 P s (k,z) , (3) 

where Tb is the mean 21-cm brightness temperature, Ps 
as before is the density power spectrum, and W/ 111 is a 
window function that expresses in Fourier space a convo¬ 
lution of density with a radial distribution function W r 
in real space. They are related by a Fourier transform 
relation: 

W[ m (k, z) = W r (r, z )~~^~dr , (4) 

where we consider a (randomly chosen) fixed point P 
at a fixed z, r is the (comoving) distance from P, and 
fo W r (r',z)dr' equals the fraction of the total X-ray in¬ 
tensity at P that comes from sources up to a distance r 
away. 

Thus, in this scenario, the power spectrum can be used 
to measure the radial distribution of X-ray flux. Be¬ 
fore we test this proposition, we note that: 1) We have 
dropped some overall (k- independent) factors, since we 
only test the shape and not the normalization of the 
power spectrum; 2) Even within the linear theory, we 
have neglected small additional complicating factors; and 
3) Non-linearity may complicate things further. 

Within our simulated cosmic volumes we measure the 
radial distribution of X-ray flux, W r . seen on average 
(i.e ., averaged over many central pixels at each redshift) 
[32l ]. How this quantity varies with radial distance de¬ 
pends on whether the radiative sources emit hard X-rays 
(in which case most of the photons travel far from the 
sources) or soft X-rays (in which case most of the pho¬ 
tons are absorbed close to the sources). We compare 
its transform IU/ m from eq. (J4J) to the following quantity 
measured from our simulated results as in mock 21-cm 
observations: 


Wk(k,z) 


PsT% 


( 5 ) 


i.e., this is the ratio of the measured isotropically- 
averaged 21-cm power spectrum to the theoretically cal¬ 


culated density power spectrum Ps (which is well deter¬ 
mined based on other cosmological probes such as the 
CMB). We have included the factors of Tj, (which is not 
directly observable) simply to make Wk dimensionless, 
but as mentioned, we are interested only in the shape of 
Wk as a function of k, and not in its overall normaliza¬ 
tion. 

Applying this method to heating fluctuations, we find 
that it is indeed possible to reconstruct the radial dis¬ 
tribution of the X-ray radiative background from 21- 
cm data (Figure [2). The slope of Wk matches that of 
Wl in remarkably well near the heating transition, i.e., 
at 2 ~ 12 — 14 for the hard SED and z ~ 14 — 18 for 
the soft SED. This works only on large scales (up to 
k = 0.1 — 0.5 Mpc -1 depending on z and the SED), 
where the heating fluctuations dominate. As expected, 
in the case of the soft X-rays the slope is flatter, i.e., 
the contribution from short distances (high k) is larger. 
The difference between the two SEDs is large and easily 
seen. Thus, measuring the slope of Wk near the heating 
transition reveals the hardness of the X-ray SED. 

Conclusions. In this Letter we have quantified for 
the first time the anisotropy in the 21-cm power spectrum 
over the wide range of cosmic epochs from cosmic dawn 
to the end of the epoch of reionization (z = 30 — 7) using 
mock 21-cm observations. We show that the anisotropy is 
large and it can be used as a cosmic “clock” to determine 
the timing of cosmic events. We also show that, in the 
case of X-ray heating, it is possible to reconstruct the 
shape of the radial flux distribution from measurements 
of the 21-cm power spectrum. This provides a new tool 
to constrain the nature of the first sources of heat in the 
Universe. 

Our findings are especially timely due to the great cur¬ 
rent interest in observations of the 21-cm signal, includ¬ 
ing the recent first constraints on cosmic heating from 
the PAPER experiment j33|. Our results will be useful 
for the interpretation of future data from radio telescopes 
such as HERA Q and the Square Kilometre Array [H], 
and should make the anisotropy of the 21-cm power spec¬ 
trum a topic of major theoretical and observational in¬ 
terest. 
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FIG. 2. A model-independent method for determining the spectrum of the X-ray sources that heated the universe. We compare 
Wk as measured from mock observations to the predicted from the radial distribution of X-ray flux as measured in our 
simulations. We compare our two cases of a hard X-ray SED (left panel) or a soft SED (right panel). We show Wk at a range 
of redshifts, z = 12 (solid red), 13 (dashed red), and 14 (dotted red) for the hard SED, and z = 14 (solid red), 16 (dashed red), 
and 18 (dotted red) for the soft SED. Each case is compared to Wl ln (which changes little with z ) at the central redshift (black 
curve); for contrast we also show Wl ln from the other panel (grey, dot-dashed line). We only wish to compare shapes, so all 
the curves are normalized to unity at k = 0.1 Mpc -1 (but to illustrate the overall amplitude at these redshifts, we note that 
the highest fluctuation at k = 0.2 is 3.3 mK for the hard SED at z — 14, and 12.4 mK for the soft SED at z — 16). We note 
that the slopes of Wk and Wk agree only on scales on which the 21-cm fluctuations are dominated by heating fluctuations; 
also note that in this Figure we are able to go to much larger scales than in Figure [l] since the isotropically-averaged power 
spectrum can be determined more accurately in a given simulation box than parameters that depend on the angular variation. 
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